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Earthquake data from Central California over the period 
January 1, 1969, to December 31, 1971, and from Western 
Nevada over the period August 31, 1954, to December 31, 1961, 
have been analyzed to investigate the hypothesis that the 
Earth tides may be the triggering mechanism for the occurrence 
of earthquakes. One method used compared the tidal components 
and the component rates -of -change at the time and location of 
actual earthquakes to those for random events uniformly 
generated over the same time period. A second method com- 
pared the differences in the Earth tides computed across the 
breadth of a fault zone at the time of actual earthquakes to 
the differences at the times of random events uniformly 
generated over the same time period. Statistical tests were 
performed to see if the distributions from the actual events 
were the same as the corresponding distributions from the 
randomly generated events. 
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I. INTRODUCTION 



Every point on the surface of the Earth is subject to 
three types of forces : the force of gravity due to the 

attraction of the mass of the Earth, the centrifugal force 
due to the rotation of the Earth, and a force due to the 
attraction of extra-terrestrial bodies, mainly the Sun and 
the Moon. The resultant forces on the surface of the Earth 
are generally referred to as the Earth tide forces. 

Because of the stresses that these Earth tides create in 
the crust of the Earth, it has been long felt that they may 
be a significant factor in the triggering of earthquakes. 
Several previous studies have been undertaken in which 
attempts were made to correlate the gravitational component 
of the Earth tides with the occurrence of earthquakes. In 
one study, Knopoff [Ref. 1] concluded that there was no 
correlation between this tidal component and the occurrence 
of the 3800 earthquakes in Southern California that he inves- 
tigated. Shlien [Ref. 2] came to a similar conclusion in his 
study of several widely separated regions of the world. 

A team of Russian researchers [Ref. 3] contend that these 
and other similar studies have suffered from two serious 
drawbacks. First, they considered only the gravitational 
component of the Earth tides. More important may be the 
tangential components which depend very strongly on the 
orientation of a fault in a region. In general, the extremes 
in the tangential components do not occur at the times of the 
extremes of the gravitational component. Second, they tried 
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to cover too large an area. Knopoff considered a region of 
18 square degrees, and the regions in Shlien's study varied 
between 50 and 450 square degrees. Since tides are time 
dependent and periodic in nature, considerations of large 
regions have averaging effects on their variation, possibly 
masking dependencies present. 
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II. OBJECTIVE 



The objective of this study was three-fold. Given a 
region of investigation and a specific time period: (1) 

compare the tidal forces from a series of actual earthquakes 
to these forces from a series of random events generated 
uniformly over the same time period, (2) compare the rates- 
of-change of the tidal forces from the actual events to those 
from the randomly generated events, and (3) compare the 
differences in each of the tidal forces over the breadth of 
the region from the actual events to those from the randomly 
generated events. These comparisons were used to test the 
hypothesis that the distribution of a force, rate, or differ- 
ence from the actual events was the same as the corresponding 
distributions from the randomly generated events. Rejection 
of the hypothesis would indicate a difference in the distri- 
butions being compared. Thus, an evaluation can be made as 
to whether a significant dependency exists between a force, 
rate, or difference and the occurrence of earthquakes for the 
particular region being studied. 
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III. ORGANIZATION 



For the purpose of this analysis, the tidal forces were 
resolved into five quantities. The orientation and con- 
struction of these quantities vvfere as follows: The gravita- 

tional component was formed along the radius of the Earth at 
the earthquake epicenter, the normal component was formed 
perpendicular to both the gravitational component and the 
strike of a fault, the axial component was formed perpendic- 
ular to the gravitational component and parallel to the 
strike of a fault, the tangent magnitude was the resultant 
formed by the combination of the magnitudes of the normal 
and axial components, and the total magnitude was the 
resultant formed by the combination of the magnitudes of the 
gravitational, normal, and axial components. By observing 
how these five quantities changed over a small change in • 
time, the rate-of -change of each quantity also was approxi- 
mated. By calculating these tidal quantities on both sides 
of a fault zone, the differences in these quantities also 
were calculated. 

It should be noted here that the foregoing definitions of 
the normal and axial components of the tidal forces usually 
precludes their being exactly tangent to the surface of the 
Earth owing to the fact that the Earth in not a perfect 
sphere. Since this angular deviation from tangency was so 
small and was nearly constant over the small range in lati- 
tude of the regions investigated, no correction was made in 
calculating component values. This correction should cancel 
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out during the comparisons of component frequency 
distributions and, therefore, have little influence on the 
results of the comparisons. 

In selecting regions for investigation, zones of high 
seismic activity were desired so that relatively small regions 
could be selected and still yield a sufficient number of 
earthquake occurrences to support a detailed analysis. The 
regions chosen were small sections of the San Andreas fault 
in Central California and the Dixie Valley fault in V/estern 
Nevada, both of which are located in the Western United 
States (Figure 1) . These two regions were chosen because of 
their differing fault structure. The San Andreas fault is a 
strike-slip fault where the land mass to the southwest side 
of the fault moves northward in relation to the land mass on 
the northeast side of the fault (Figure 2) . In the partic- 
ular region of the San Andreas fault investigated, the strike 
of the fault is approximately 40 degrees west of north as 
measured from longitude 121 W. The Dixie Valley fault is a 
graben fault where the material in the center of the fault 
zone moves downward in relation to the sides of the fault 
(Figure 3) . The strike of this fault is approximately 
parallel to longitude 118 W. Co-ordinates of these regions 
are given in Table 1. 

In choosing a time period over which earthquakes were to 
be analyzed, a large enough time period was needed so that 
bias due to the periodic fluctuation of the Earth tides 
(monthly, annually, etc.) would be damped or averaged out. 
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This was needed so that differences found in the comparison 
of these quantities xvould not be attributable to the cyclic 
nature of the tides. For the specified portion of the San 
Andreas fault, 1732 earthquakes of Richter magnitude 0.0 or 
greater occurred in the three year period from January 1, 

1969, to December 31, 1971. For the portion of the Dixie 
Valley fault, the available data prior to 1962 was for shocks 
of magnitude 4.0 or greater, and from 1962 on, for shocks of 
magnitude 0.0 or greater. So as not to unduly bias the 
distribution of data with the increased frequency of recorded 
small shocks after 1962, only the time period through 1961 
was considered. During the seven year period from August 31, 
1954 (the start of detailed recording in this area) , to 
December 31, 1961, 100 earthquakes of magnitude 4.0 or greater 
occurred . 
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Figure 1 
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San Andreas 

Strike-Slip Fault Structure 
Figure 2 




Dixie Valley 
Graben Fault Structure 

Figure 3 
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IV. EARTH TIDE DISCUSSION 



The actual tidal data were generated through the use of a 
computer routine (unpublished) by Professor R. H. Shudde, 

Naval Postgraduate School. Much work has been accomplished 
over the past few years on the subject of theoretical Earth 
tides. Melchior [Ref. 4] has one of the better discussions 
on the subject. He goes into great detail concerning the 
computation of the theoretical Earth tides using an equi- 
potential surface argument. His procedure determines the 
locations of the Sun and the Moon in relation to the Earth at 
a particular time and establishes the potential surfaces for 
the tides. Then by specifying the location of a point on the 
Earth in relation to these surfaces, he determines the tides 
at that location. Shudde 's program takes a different approach. 
It also determines the locations of the Sun and the Moon using 
the procedures in reference 5, but then at a particular point 
on the Earth, the tides are calculated directly using the 
standard gravitational force laws. For a more detailed 
description of the Earth tide calculations, see the discussion 
in Appendix A. It will suffice here to say that all Shudde 's 
program requires for inputs are the time, date, and geographic 
co-ordinates of an earthquake occurrence. For this analysis, 
the times and locations for the actual occurrences were taken 
from references 6-8 and from data supplied by the National 
Center for Earthquake Research. The dates for the random 
occurrences \vere generated uniformly over the time periods 
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specified in Section III. The location for these random 
events was taken to be the center of the region being 
investigated . 
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V. ANALYSIS 



The objective of this analysis was to test the hypothesis, 

H^: A = R 

where A is the distribution of data from the actual events, 
and R is the distribution of the corresponding data from the 
randomly generated events. Rejection of this hypothesis for 
a particular force, rate, or difference would indicate a 
difference in the two distributions being compared and, hence, 
a dependency between that quantity and the occurrence of 
earthquakes . 

The primary statistical test used to test this hypothesis 

was the Kolmogorov-Smirnov (K-S) , Two-Sample test. This is a 

non-parametric test independent of distributional assumptions 

about the data other than the distribution be continuous. It 

is an exact test and operates by the comparison with the 

actual data rather than with grouped data. A second test 

used was the Chi-Square Goodness -of -Fit test which operates 

by the comparison of grouped data. 

The K-S statistic "possesses the advantage that it is an 
exact method, whereas the Chi-Square method requires a 
fairly large sample to justify the approximations that are 
needed to apply it" [Hoel, 1971]. 

The Chi-Square test was used because this is the common 

statistical test used in many studies. The K-S test was used 

because it is a more powerful test and was more appropriate 

to the data. 

A. KOLMOGOROV-SMIRNOV, TWO -SAMPLE TEST 

This test was implemented through the use of the IBM 
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computer subroutine, K0LM2 [Ref. 10], with a modified sort 
routine [Ref. 11]. Using these routines, for each of the 
five tidal, rate, and difference quantities, the distribution 
of data from the actual earthquakes in a region was compared 
to the distribution for that quantity calculated from the 
randomly generated events. The two sets of data for each 
distribution were each sorted into increasing order and formed 
into sample distribution functions. The maximum absolute 
ordinate difference between the two distributions was deter- 
mined, A sample size adjustment factor was calculated by 
taking the square root of the ratio of the sample size 
product to the sample size sum. The K-S, two-sample statistic, 
Z, was then calculated by taking the product of the maximum 
difference and the adjustment factor. That is. 



where A(x) is the distribution function from the actual events, 
S(x) is the distribution function from the random events, 
is the number of events in A(x), and M 2 is the number of 
events in S(x). The use of the adjustment factor is necessary 
in the K-S test. If not applied to the statistic as shown 
above, the inverse of the factor would need to be applied to 
the critical values listed in Table 2 to keep the test statis- 
tic consistent with the critical values. In this case, appli- 
cation of the adjustment factor to the statistic simplified 
the use of the computer routines and enabled the critical 
values to remain constant for a given level of significance. 
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B. CHI-SQUARE GOODNESS-OF-FIT TEST 

For each tidal quantity, the maximum range covering both 
the tidal data from the actual earthquakes and the data from 
the randomly generated events were divided into 20 equal 
intervals and the data was grouped into these intervals . 
Adjacent intervals \^ere collected together until the expected 
frequency in the collected intervals exceeded three [Ostle, 
1963] . The Chi-Square statistic was then calculated as 




where 0^^ is the frequency in interval, i, from the actual 
events, E^j^ is the expected frequency in interval, i, and m' 
is the number of collected intervals. Selected critical 
values for the Chi-Square test are given in Table 3. 
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Level of 
Significance 


Value of Z So Large as to 
Call for Rejection of the 


.10 


1.22 


.05 


1.36 


.025 


1.48 


.01 


1.63 


.005 


1.73 



Adapted from Beyer, p. 429. 

Critical Values of the 
K-S, Two-Sample, Statistic, Z 

Table 2 
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Adapted from Beyer, p. 294. 



Critical Values of the 
Chi-Square Statistic 

Table 3 
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VI . RI'SULTS 



In the analysis of the tidal data, some bias in the 
distributions of the data was expected due to the physical 
distance between the origin of the actual earthquakes and the 
geometric center of the region at which the random events 
were generated. Although this bias was felt to be minimal 
due to the small size of the regions and the large number of 
events, additional checks were performed to test the size of 
this bias. Random events were generated at the north, south, 
west, and east extremes of each region and the five tidal 
and five rate quantities were calculated at each location for 
each occurrence. The ten quantities at the north and south 
extremes were then compared using the K-S test, as were the 
ten quantities from the west and east extremes. The results 
are shown in Tables 4 and 5. The conclusion was that bias in 
the data due to the physical separation of the events was so 
slight as to have an insignificant influence on the final 
analysis. It must be stressed here that the significance of 
the following results pertains only to the particular region 
and the nature and orientation of the particular fault being 
analyzed . 

A. SAN ANDREAS 

The K-S, two-sample statistics for each of the five 
tidal, five rate, and five difference quantities acting in 
the San Andreas region are shown in Table 6 and the Chi- 
Square statistics for these quantities are shown in Table 7. 
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Referring to Tables 2 and 3 for a 0.05 level of significance, 
the normal tidal and gravitational rate quantities were 
statistically significant under both tests. None of the 
tidal differences was statistically significant at this level 
in either test. 

B. DIXIE VALLEY 

The K-S, two-sample statistics for each of the five tidal, 
five rate, and five difference quantities acting in the Dixie 
Valley region are shown in Table 8 and the Chi-Square statis- 
tics for those quantities are shown in Table 9. Referring to 
Tables 2 and 3 for a 0.05 level of significance, the normal 
rate quantity was statistically significant under both tests. 
None of the tidal difference quantities was significant at 
this level in either test. 

After much of this analysis had been completed, a 
similarity between the methods of calculating the rate and 
difference quantities was observed. The rate quantities were 
calculated by computing the tidal quantities at a location, 
allowing a small unit of time to elapse, and then computing 
the tidal quantities again. The change in the tidal quantity 
divided by the change in time was used as its rate of change. 
Because of the rotation of the Earth, changes in time corres- 
pond to changes in location at a fixed latitude. A lapse of 
four minutes is equivalent to a distance change of one degree 
in longitude. For small time lapses, the distance considered 
moved is negligible. For larger time lapses, the calculations 
can be considered as being at two separated locations. The 
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difference quantities were calculated by computing simultan- 
eously the tidal quantities at tv'^o locations separated by one 
degree in longitude. Since both procedures involve the 
calculation of components at separated locations, there 
should have been a correlation bet\\^een the test results on 
the rate quantities and the test tesults on the difference 
quantities. Comparison of the San Andreas results in Table 
6 showed some similarity but the interpretation was clouded 
in that the tangential components in this region were rotated 
through an angle of 40 degrees in order to align them with 
the strike of the fault. For the region in which the 
tangential components had not been rotated, comparison of the 
Dixie Valley results in Table 8 showed a marked agreement in 
the order of the significance of the results. Therefore, 
these interpretations of the results could have been made 
from a single analysis of the rate quantities without doing 
the difference calculations. 
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Component 



Component 


K-S Statistic 
N-S W-E 


Gravitational Tide 


.77 .44 


Normal Tide 


.33 .25 


Axial Tide 


.39 .28 


Tangent Magnitude 


.32 .23 


Total Magnitude 


.39 .25 


Gravitational Tide Rate 


.42 .36 


Normal Tide Rate 


.21 .20 


Axial Tide Rate 


.25 .23 


Tangent Magnitude Rate 


.34 .29 


Total Magnitude Rate 


.53 .42 



Results of Bias Check in Tidal Comparisons 
N-S and W-E Extremes of the San Andreas Area 
Using the K-S Statistic 

Table 4 
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Component 



K-S Statistic 
N-S W-E 



Gravitational Tide .91 .17 



Normal Tide 




.28 


. 16 


Axial Tide 




.48 


.20 


Tangent Magnitude 




.32 


.16 


Total Magnitude 




.36 


.15 


Gravitational Tide Rate 




.46 


.17 


Normal Tide Rate 




.29 


.17 


Axial Tide Rate 




.26 


.15 


Tangent Magnitude Rate 




.30 


.20 


Total Magnitude Rate 




.61 


.16 



Results of Bias Check in Tidal Comparisons 
N-S and W-E Extremes of the Dixie Valley Area 
Using the K-S Statistic 

Table 5 
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Component 


Statistic 


Gravitational Tide 


1.1411 


Normal Tide 


1.5038 


Axial Tide 


1.1776 


Tangent Magnitude 


1.3344 


Total Magnitude 


0.8658 


Gravitational Tide Rate 


1.6515 


Normal Tide Rate 


0.6126 


Axial Tide Rate 


0.8469 


Tangent Magnitude Rate 


0.9541 


Total Magnitude Rate 


1.2273 


Gravitational Tide Difference 


1.1724 


Normal Tide Difference 


1.2752 


Axial Tide Difference 


0.9117 


Tangent Magnitude Difference 


1.2259 


Total Magnitude Difference 


1.1985 



Results of Actual Tidal 
San Andreas Area Using 


Comparisons in the 
the K-S Statistic 



Table 6 
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Component 


Statistic 


Degrees o£ 


Gravitational Tide 


40.22 


19 


Normal Tide 


96.84 


18 


Axial Tide 


36.99 


17 


Tangent Magnitude 


23.37 


18 


Total Magnitude 


34.33 


18 


Gravitational Tide Rate 


114.21 


19 


Normal Tide Rate 


38.38 


18 


Axial Tide Rate 


23.48 


17 


Tangent Magnitude Rate 


38.21 


18 


Total Magnitude Rate 


81.64 


17 



Results o£ Actual Tidal Comparisons in the 
San Andreas Area Using the Chi-Square Statistic 

Table 7 
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, 



Component 



Statistic 



Gravitational Tide 


1.2218 


Normal Tide 


0.6278 


Axial Tide 


0.8792 


Tangent Magnitude 


0.7129 


Total Magnitude 


0.6139 


Gravitational Tide Rate 


0.6535 


Normal Tide Rate 


1.4119 


Axial Tide Rate 


0.5901 


Tangent Magnitude Rate 


0.4673 


Total Magnitude Rate 


0.9228 


Gravitational Tide Difference 


0.5743 


Normal Tide Difference 


1.2654 


Axial Tide Difference 


0.6040 


Tangent Magnitude Difference 


0.5822 


Total Magnitude Difference 


1.1605 



Results of Actual Tidal Comparisons in the 
Dixie Valley Area Using the K-S Statistic 

Table 8 
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Degrees of Freedom 



Component 

Gravitational Tide 
Normal Tide 
Axial Tide 
Tangent Magnitude 
Total Magnitude 

Gravitational Tide Rate 
Normal Tide Rate 
Axial Tide Rate 
Tangent Magnitude Rate 
Total Magnitude Rate 



Statistic 



11.22 13 

6.90 14 

26.30 14 

15.04 15 

10.47 11 

7.34 14 

25.47 14 

5.99 15 

3.91 12 

13.72 13 



Results of Actual Tidal Comparisons in the 
Dixie Valley Area Using the Chi-Square Statistic 

Table 9 
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VII. CONCLUSIONS 



Results of the analysis indicate that, even at high levels 
of significance, there exist components of the Earth tides 
and their rates -of -change that are significantly related to 
the occurrence of earthquakes and that these components vary 
with the nature and orientation of the fault being investi- 
gated. 

The significance of the normal tidal component on the 
strike-slip San Andreas fault gives credence to the theory 
that the occurrence of earthquakes in this fault zone may 
result from the lessening of the frictional force acting 
between the two. vertical faces of the fault. Since the 
frictional force acting between two surfaces is a function of 
the normal force acting on those surfaces, a lessening of the 
normal force could decrease the frictional force enough to 
allow the accumulated stress in the fault to be released. 

The relationship bet\\reen the rate -of - change of a tidal 
force and the difference in a force over the breadth of a 
fault has been discussed previously. Although no difference 
quantities were considered significant at a 0.05 level of 
significance, reference to Tables 2, 6 and 8 reveal that these 
quantities did differ considerably and could have been consid- 
ered significant at loiver levels of significance. The signif- 
icance of the gravitational rate in the region of the San 
Andreas fault indicates a differential in the vertical forces 
acting on opposite sides of the fault. This significance is 
difficult to interpret in view of the strike-slip nature of 
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the fault. Because of the proximity of this fault to the 
Pacific Ocean, this significance may indicate the presence of 
ocean tidal loading effects on the occurrence of earthquakes 
in this region. The significance of the normal rate in the 
region of the Dixie Valley fault indicates a differential in 
the east-west force acting on opposite sides of this fault. 
This interpretation is consistent v^?ith the nature of the 
Dixie Valley fault in that the downv\^ard slippage of the center 
of the fault zone occurs with the separation of the sides 
of the fault. 
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APPENDIX A 



A. OBSERVATIONS ON ASSOCIATED LEGENDRE POLYNOMIALS. 



The generating function for the associated Legendre 
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B. EARTH TIDE CALCULATIONS, 




R = Earth radius vector to epicenter 
a = Scalar magnitude of R 

R- = Distance vector from epicenter to B- 
r- = Scalar magnitude of R- 

C ““J* c 

R^q = Distance vector from origin to B^ 
r^Q = Scalar magnitude of R^^^ 

B^ = Extra-terrestrial body, i=l (moon), 2 (Sun) 
F- = Force vector between the origin and B- 
F- = Force vector between the origin and B- 
= Angle between R and Rj^^ 
m^ = mass of B^ 



From Newton’s laws of gravitation. 



Fie- 



u. ^ 



C o 



V- ^ 

'i, o 



The total tidal force vector, T, is the sum over all bodies 



T .1 

(La I 



where Tj^ = Fj^Q - Fj^g = k*mj^ 



^io ^ie 



r . ^ r • ^ 
10 le 



Since 



k* rv,,: 









T' -~r 

With kje = + <3-'’ - Z a.' ir^o cosC°^t), the expressi 

be written as 



V"i. 



on 






(7) 



can 
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2(— |Cos6<c) 




3^ 

Z 



(?) 



Letting x^ = cos (o<-L ) , and y^ = 



10 



, a comparison of (4) with 



-i" 



(8) yields, 

V^'ic/ 5=' 

^ Fc'(cosC<^^)) 

= ^ Ir.o/ 

c^. 



5 t M 






-1 a) 

9'-' 









d cosCot,\^ 



Substituting (9) into (7) yields, 

o , 4 p ('cosfa.')) 

Rio -- (Eoo -R) Z (lr,oj j cos(!o<;) * 

fi=' 



Ti = 



k ' Wit 



to L 



Letting tj^ be the component of in the direction of N, a 
unit vector in the plane normal to R, and letting be the 

i# 

angle between N and , results in the general expression, 
^n~— *— i’ after rearranging terms. 



k.«; fn r 'i'' “* 

^ t,o ^ 



cosC^i) 



(co5C'=<i)) 0®) 



Expanding terms in (10) and substituting values from (5) , 



K • 



~l — Ctss 



t,o 



(ftl) 






4 



<1 co^(o<:) 



Tz (^C os Ce< 






4^'®' <4 cosCoii) 

- ^ ccs(/3i,)j^3co5Co<0+ j)**""* (") 

R‘R- X. 

Noting 'that cos(o<i)= — = , 

a • r ■ r - 
10 10 

if N is taken to be the north(z) direction, (11) gives t 
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t 



w 



ith cos(|3;) = 



z ; 



lO 



, and i£ N is taken to be the east(y) 

Yi 



direction, (11) gives t with cos(/^l)= — If the 

^ ^io 



co-ordinate system is oriented so that R is in the x-direc- 

R-T. 

tion, then t-^-=t = 



— —1 



— ^ ^ c 



cx*)r. 



.3 



CO 



a • - Co~‘ -cosC^O -<^) 

<J cotfoO Pa(.'^°sC"<0)J j 



<» , a \ 

+ I (-J 






I 

L 



k • vv); 






i-f 



Expanding a few terms of (12) , collecting coefficients of 
j-^y ^ where j=0,l, and substituting values from (6), 



with cos(c<i)= - — • For tangential components in other than 
^io 

the y and z directions, N can be any unit vector in the 
plane normal to the vector R. 
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